Probing quantum coherence in qubit arrays 
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We discuss how the observation of population localization effects in periodically driven systems 
can be used to quantify the presence of quantum coherence in interacting qubit arrays. Essential 
for our proposal is the fact that these localization effects persist beyond tight-binding Hamiltonian 
models. This result is of special practical relevance in those situations where direct system probing 
using tomographic schemes becomes infeasible beyond a very small number of qubits. As a proof 
of principle, we analyze a scheme consisting of a chain of superconducting flux qubits and show 
that localization effects underpinned by quantum coherent interactions should be observable within 
realistic parameter regimes. 



I. INTRODUCTION 

Transport processes are of fundamental importance in 
a wide variety of physical and biological systems, ranging 
from the actual motion of particles on a lattice [TJ [2] to 
the transfer of classical and quantum information across 
spin or harmonic chains [3H5]. Relevant for our pur- 
poses, there exist specific features of the transport pro- 
cess that are intrinsically linked to the dynamics of the 
chain and in particular to whether or not the chain el- 
ements can interact coherently [B]. In the mid-80's the 
motion of a charged particle on a onc-dimcnsional lat- 
tice under the influence of a time-dependent electric field 
was studied and shown to exhibit dynamic localization 
(DL) pp. The canonical situation to illustrate this phe- 
nomenon is provided by an infinite linear chain of sites 
along which a charged particle moves under the combined 
influence of a nearest-neighbour exchange interaction and 
a time-dependent external driving. In that setting, it was 
found that the mean-square displacement of the particle 
as a function of the field modulation E\, rather than 
exhibiting a diffusive behaviour, does not grow without 
bounds but oscillates sinusoidally. A related phenomenon 
is the so called coherent destruction of tunneling (CDT), 
initially formulated for a symmetric, externally driven, 
double- well potential [7]. Both DL and CDT are gen- 
uine manifestations of coherent quantum effects resulting 
from the interference between different transition paths 
that leads to the selective inhibition of transport j8]. In 
contrast, in the classical case, and from an initially lo- 
calized state, an equilibrium state would be attained in 
which neighbouring sites would be equally populated. 
Ample experimental evidence supports the existence of 
both types of localization effects in a variety of systems. 
DL has been observed in Rydberg atoms, where the lo- 
calization regime is characterized by a " freezing" of the 
width of the wave packet with respect to the Rydberg lev- 
els [5]; in driven quantum wells or semiconductor super- 



lattices, where a suppression of the conductance was ob- 
served [10] ; in ultracold atoms interacting with a standing 
wave of near-resonant light, where this phenomenon was 
found in the suppression of momentum jTTJ [T^] . CDT 
has also been recently observed in both noninteracting 
[T3] and interacting systems [13] . Further, motivated by 
the desire to study the effects of quantum coherence and 
dephasing noise and the interplay of the two on transport 
processes in biological systems [TS], recently the detec- 
tion of dynamic localization was proposed as a way to 
demonstrate the possible existence of coherence effects 
in ion channels (16] , i.e. protein complexes that regulate 
the flow of particular ions across the cell membrane and 
that are essential for a wide variety of cellular functions. 

Here we extend this work and discuss the possibil- 
ity of observing localization effects beyond the canonical 
setting, including deviations from a strict tight-binding 
Hamiltonian as well as the inclusion of non-Hamiltonian 
(noisy) effects. We will show that signatures of localiza- 
tion effects can still be observed in this case and apply 
these results to the problem of qualitatively witnessing 
quantum coherence in an interacting chain of supercon- 
ducting qubits. 

II. RENORMALIZATION OF INTRA-QUBIT 
INTERACTIONS BY MEANS OF AN 
EXTERNAL MODULATION 

Motivated by specific qubit realizations in the solid 
state, we analyze an array of interacting qubits subject 
to a Hamiltonian of the form {K = 1): 

N N-l 
k=l k=l 

with N denoting the number of qubits in the chain, 
the site energies for each qubit, and Jk,k+i the coupling 
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between neighbouring qubits k and k + 1. In the presence 
of a time-dependent external driving of the form 



1 N 

H ac (t) = - 2J k ■ E ac cos(ut)al 



(2) 



k=l 



the Hamiltonian i/(i) = Hq + H ac (t) of the chain reads 
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(3) 

With the substitution crjj; = cr^T + a k , the Hamiltonian 
above can be rewritten as the sum of three contributions, 



with, 



H(t) = H z (t) + H 1 +H 2 
1 N 

H z (t) = - Y (^k + kE ac cos(ujt))a z k 
fc=i 
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H 2 = Y Jk,k+i{(ri^k+i +h.c). 
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Dchning the total excitation number operator as 

JV 



(8) 
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it is easy to see that [H z {t),N] = [i/i,iV] = while 
[H.2,N\ 0. For this reason, the term H\ is usually 



referred to as an exchange interaction, in the sense 
that it allows for a hopping of the excitations within 
the chain, but it does not create nor annihilate them. 
This is the canonical interaction in previous studies of 
dynamical localization in systems that can be modelled 
with a tight-binding Hamiltonian pQ. 

In the following lines we will show that the interac- 
tions described by the terms H\ and H2 can be indeed 
enhanced or inhibited separately by the proper tuning of 
frequency lo and amplitude E ac of the external field. 
To gain insight into the problem, it is convenient to move 
to an interaction picture with respect the time-dependent 
term H z (t). That is, we first define 



U (t) ee exp + |~*y rfT_ff 2 (r)| = exp j-i J <1tH z (t) 

(9)' 

where we have made use that H z (t) commutes with it- 
self at different times to write the last equality above. 
Computing explicitly the integral above and taking into 
account that the operators acting on different sites com- 
mute, we have that 
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ui k t .kE ac . 



(10) 



Hence, the interaction picture Hamiltonian of our chain 
can be written as 

H\t) ee U Q {t)\H t + H 2 )U Q (t) ee H[(t) + H' 2 (t), (11) 

where we have defined H[(t) = Uo(tyHiUo(t) and 
H'Jt) ee Uo(t)^H 2 Uo(t). Now, making use of eq. 
(10 1 and the Jacobi- Anger expansion e tzsln (4>) — 



12n=-oo Jn{z) e%n i where J n (x) are the Bessel functions 
of first kind, we can proceed further and evaluate the 
form of these terms explicitly as: 
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For the sake of clarity we will consider in the follow- ing the case of a homogeneous chain with 0Jk = ujq and 
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Jk,k+i = J f° r a ll values of fc. Under this assumption the 
Hamiltonian H'(t) can be written as: 



JV-l 



H'(t) = ^ .g(t)<7+ ffH1 + g' k (t)a+a+ +1 + h. c. (14) 



k=l 



with time-dependent rcnormalized couplings g(t) and 
g'(i) defined as 



g(t) = J ]T J n 



n— — oo 



(15) 



(16) 



In the regime where the tunneling frequency of the qubits 
is much smaller than the frequency of the driving field, 
that is J«w, we can invoke the rotating wave approxi- 
mation in the series above and neglect those terms that 
rotate faster than J. In particular, for g(t) this means 
that only the non-rotating term with n = survives and 
we can write 



controlled manner and in a variety of geometries makes 
them a promising candidate for the implementation of 
quantum registers and information processers |17) . On 
the other hand, while the fabrication of structures in- 
volving many qubits is indeed feasible [TS], its effective 
probing, and in particular the verification that the sys- 
tem does exhibit quantum coherence, is beyond the realm 
of current technology even for moderate system size. The 
current state of the art is provided by the tomographic 
analysis and the entanglement verification of 3 qubit sys- 
tems HH2T]. 
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FIG. 1: Schematic depiction of a superconducting qubit chain 
with controllable injection and readout as well as individ- 
ual addressing for local modulation. By means of population 
measurements on a selected qubit, the presence of quantum 
coherence is inferred from dynamical localization effects, as 
described in the main text. 



i(t) = g = J ■ Jo 



(17) 



Applying the same reasoning to g'(t) it follows that the 
only possibility to have surviving terms is that a reso- 
nance between ujq and u ocurrs such that 2wo + n'u) = 
for some integer value n! G Z. In this case, eq.(16) can 
be further simplified, yielding 



j{t) = g' = J- J M 



\n'\E a 
2u 



(18) 



We therefore see that the effect of the external modu- 
lation can be interpreted as renormalization of the cou- 
pling constants, imprinting a periodic dependence that 
will lead to a selective inhibition of transport. In the fol- 
lowing sections we will provide numerical evidence of the 
accuracy of the expressions derived above. This type of 
localization effect will be later exploited to detect signa- 
tures of coherent interaction in arrays of coupled qubits. 



III. SYSTEM DESCRIPTION 

We will study the persistence of population localiza- 
tion effects induced by the renormalization of the hop- 
ping coupling, as explained in the previous section, in 
the dynamics of superconducting qubit arrays. Supercon- 
ducting qubits are effective two-level systems with a con- 
trollable transition frequency between their eigenstates, 
whose potential to be manufactured lithographically in a 



In general the superconducting qubit Hamiltonian can 
be written as H s = (e/2)<j z + (A/2)a x . Depending on 
the particular qubit realization the parameters e and A 
refer to different variables defining "charge", "phase" or 
"flux" qubits. The latter, also called "persistent current 
qubit" , consists of a superconducting loop interrupted 
by three Josephson junctions, two with capacitance C\ 
and the third with C 2 [H (221 E3]. The values of the 
three Josephson junctions coupling constants, Ej\ cor- 
responding to capacitance C\ and Ej^ to capacitance 
Ci respectively, are chosen so that the Josephson part of 
the Hamiltonian alone defines a bistable system. At a 
value of external magnetic flux $ m 0.5<E>o (^o = h/2e 
is the superconducting flux quantum) , the system carries 
either a clockwise or counterclockwise persistent current, 
each generating an equal but opposite magnetic flux and 
defining the two possible states of this qubit. With an 
appropriate choice of the parameters -Ej,i,2 and C\$, the 
barrier in phase space separating the left and right cur- 
rent states can be made low enough so that tunneling 
between the two classical states can occur. For the flux 
qubits A represents the tunneling amplitude, while the 
energy bias e = — ip(3> — 0.5<I>o) is proportional to the 
detuning $ — 0.5$oi with I p the circulating current. The 
energies of the ground and first excited states are thus: 



E^ 



+ A 2 . 



For flux qubits, the most natural implementation for 
the interaction is through the mutual inductances be- 
tween the loops. The flux generated by one qubit loop, 
which depends on its internal state, adds to the total flux 
picked up by the neighbouring qubits, thereby changing 
the energy biases of those qubits. The mutual indue- 
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tance, and therefore the strength of the interaction (J), 
depends on the geometry of the qubit loops, specifically 
their size and their proximity to each other. The strength 
of the coupling can be enhanced in two steps. First by 
physically connecting the loops so that their persistent 
currents share a common line. In this case the kinetic 
inductance of the shared line adds to the geometrical in- 
ductance, where the former term can easily be the dom- 
inating part in the total interaction strength. To reach 
very strong coupling, for example to reach the regime 
J > A, a fourth junction can be placed in the shared 
line. If the capacitance and Josephson junction coupling 
constant of this fourth junction are large compared to 
parameters of the qubit junctions, the single-qubit prop- 
erties are not significantly altered, while still the inter- 
action strength can be enhanced dramatically. Coupling 
strengths of several GHz are easily achieved. For the sake 
of the suggested experiments in this paper it is however 
not necessary to reach such high coupling values since 
the hopping inhibition between the qubits is enhanced in 
the regime where the hopping constant is smaller than 
the driving frequency i.e. J«w (see eqs. (15 16 1 and 
explanations below those expressions). 

In the experiment described in reference |24j . it is 
demonstrated how the standard flux qubit design, for 
which only e is a tunable parameter, can be extended 
to also have a tunable A. The main change is that the 
smallest junction of the qubit is replaced by a small loop 
containing one junction in each arm, i.e. in a SQUID 
(Superconducting Quantum Interference device) geome- 
try. The SQUID acts as a single junctions with tunable 
Ej, tuned by the flux penetrating the small loop. Lo- 
cal control lines can be used to change the flux through 
this loop, thereby controlling the Ej2 of the qubit, and 
thus the A. Fig. [T] shows a schematic of a possible imple- 
mentation of a chain of strongly interacting flux qubits 
with tunable tunneling splitting A. When the array of in- 
teracting superconducting flux qubits is operated at their 
corresponding degeneracy point the chain is described by 
a Hamiltonian of the form eq.([T]), with u)i = Aj. Typical 
values for the tunneling amplitudes are in the range of 5- 
20 GHz while nearest neighbour couplings J ~ 200 MHz. 
The residual next-to-nearest coupling is smaller than 10 
MHz. 

As far as the coupling to the environment is concerned, 
typically, noise sources that are located relatively far 
away from the chain couple to multiple or all qubits in 
the chain, for example magnetic coils, or some parts of 
the control and readout circuits. Noise sources that are 
located much closer, such as local control lines for the 
individual qubits, or microscopic noise sources in the ma- 
terials surrounding the qubits, couple only, or mostly, to 
a single qubit. In this work we focus on the latter type: 
we suppose the array to be in contact with an external 
environment that acts locally on each qubit and can lead 
in principle to both pure dephasing and dissipation. In 
our model each qubit is coupled to its local environment 
via a spin-boson Hamiltonian of the form 



H. 



SB 



He, + Hh + H a —h 



(19) 



where the bath is modelled as a collection of harmonic 
oscillators, Hi, — J^ fe cj)-& k o,k and system and bath couple 
linearly through the Hamiltonian [25] 



crlX, 



(20) 



Here Jtj = J2k9k(^k + %) denotes the bath's force op- 
erator. In the qubit eigenbasis and using the same sym- 
bols to denote the Pauli matrices in this new frame in 
order not to complicate notation, we can derive, using 
the standard assumptions, a Markovian master equation 
p = -i[H,p] + Ldeph(e) + L diss (g) for the qubit array 
where the coupling to the local environment is described 
in terms of Lindblad terms of the form 



N 

idoph(g) = Tdoph y2(2a+ o-rga+o-r - {afar,g}) (21) 



£diss(p) = 7d 



N 

>ssJ2(2o-rga+ 



(22) 



with 7de P h and 7di ss denoting the dephasing and dissipa- 
tion rate, respectively, af — (of ± «of )/2 acting on the 
i-th qubit of the chain and { , } the anticommutation op- 
eration. The value of the noise rates in these expressions 
depends strongly on the selected qubit operating point 
via the parameter 8 — arctan(A/e). Measured values for 
Ti = 7jJ s range from 150 to 500 ns while pure dephasing 
times (TJ ) are typically around 300 ns. In the specific 
case where each qubit in the chain is operated at the de- 
generacy point, so that e = for every qubit, the pure 
dephasing term, which has a rate 7d ph ~ cos(#) [26], 
cancels out and the chain is subject to dissipative noise 
only. This is the first parameter regime that we are going 
to analyze in the next section with the aim of unveiling 
coherent dynamics through dynamical localization effects 
in a driven chain. Later, still operating each qubit at its 
degeneracy point, we will relax the constraint of having 
only dissipation noise and we will consider possible de- 
phasing effects arising from terms with the form of eq. 
([211. 



IV. NUMERICAL RESULTS 

The phenomenon described above concerning the 
renormalization of the coupling constants in the inter- 
action Hamiltonian eq.Q under the effect of an exter- 
nal sinusoidal driving field can be clearly observed in 
fig. [2] There we study the dynamics of a chain consist- 
ing in two inductively coupled superconducting qubits 
operated at their degeneracy points, so that the (un- 
driven) system Hamiltonian is given by eq.Q. In this 
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FIG. 2: Coherent destruction of tunneling. Simulations 
corresponding to a N = 2 chain with wi = ui2 = 10GHz and 
hopping J = 10MHz. The black dots have been computed 
with to = 0.3GHz and H(t) = H z (t) + Hi (see definitions 
in the text). The black line has been computed with u = 
0.3GHz and H(t) = H z (t) + Hi + H 2 . The red line has been 
computed with lo = 2.0GHz and H(t) = H z (t) + Hi + H 2 . 
The dissipation rate is the same in all curves 7di ss = 1MHz. 
The dashed vertical lines represent the first zo,i — 2.4048 and 
second 20,2 — 5.5201 zeros of the Bessel function Jo(z). 



Case, as discussed before, the noise is purely dissipative 
when expressed in the rotated basis. The system is ini- 
tialized so that only the first site is excited. The different 
curves correspond to the maximum value of the popula- 
tion that has been transferred from the first to the sec- 
ond qubit within a sufficiently long time interval (1 /is). 
The black line corresponds to an off-resonance situation 
(ljq = 10 GHz, u> — 0.3GHz, with no integer value n' such 
that 2wo + n'uj — 0) . In this situation the coupling of the 
contribution given by H% effectively renormalizes to zero 
and the dynamics is governed by the tight-binding term 
Hi . According to expression eq. ( 15 ) , it is the Bessel func- 



tion Jo(E ac /u}) which governs this behavior and when 
its argument E ac /ui coincides with one of its zeros, then 
the hopping between both qubits is suppressed. On the 
other hand, the red curve has been computed on a res- 
onance situation (luq — 10 GHz, ui = 2.0GHz, such that 
2u)q + n'oj = for n' — — 10^ In this situation both 



couplings g and g' in eqs.(15 161 arc in general different 



from zero and the total dynamics is hence more convo- 
luted. Notice that the fact that both the red and black 
curves are indistinguishable to the eye for low values of 
E ac /uj is due to the slow buildup of the Bessel function 
J n i (z) (with n' = 10) that governs the H2 term. Finally, 
the black circles in fig. [2] have been computed using the 
same parameters that we used for the off-resonance situa- 
tion but this time we explicitly excluded the contribution 
given by the H2 terms. The fact that the black line and 
the black circles superimpose each other is a clear indi- 
cation that the term H2 effectively renormalizes to zero 
when the field is out of resonance with this term. 

The existence of resonance conditions for the terms 



FIG. 3: Resonance conditions. Simulations for a N = 2 
chain with wi = u)2 — 10GHz, J = 10MHz and a dissipation 
rate 7di ss = 1MHz. For each value of u), the electric field is 
fixed so that E &c /uj = zo,i, where zq.i — 2.4048 is the first 
zero of the Bessel function Jo(z). The axis on the top repre- 
sents some integer values n' — 2uuo for which the resonance 
condition in the coupling constant g' is fulfilled. 
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FIG. 4: Effect of frequency inhomogeneities. Maximum 
value of the population transferred from qubit 1 to qubit 2 in a 
N = 2 chain after a fixed time i ma x = 500ns. The frequency of 
the external field is ui — 1.3GHz and its magnitude is fixed so 
that E&c/uj = 1.2, far from the closest zero of the Jo(z) Bessel 
function. The parameters of the qubits are ujq = 10GHz and 
J — 10MHz. The dephasing and dissipation rates have been 
set to zero for clarity. 



contributed by H2 in expression eq.Q is clearly il- 
lustrated in fig. [3j which has been evaluated with the 
same parameters as in fig. [2j We have however fixed 
the ratio E ac /uj — zo,i with zq,i — 2.4048 the first 
zero of the Bessel function of order zero. With this 
constraint the population transfer induced by the terms 
contributed by Hi in eq.Q is completely suppressed. 
We have however the freedom to use different values 
of the driving field oj such that the term H2 is either 
on- or off- resonance with the frequency ujq. We can 
clearly see in fig. [3] those values (compatible with the 
discrete grid used to sample the frequency u>) where the 
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FIG. 5: Effect of pure dephasing noise. (top) Popula- 
tion transfer from qubit 1 to qubit 2 in a N = 2 chain and 
the following parameters: ui\ — u)2 = 10GHz, J = 10MHz, 
oj — 1.3GHz. The dissipation rate has been set to zero (7di ss = 
0.0) for clarity and a new dephasing rate jdeph has been in- 
cluded (£ d eph = IdcphS-Cfc^ 2 ), with £ k (a z ) = (J'tpal - p). 
Each panel corresponds to a fixed value z = E ac /^ with 
2o,i — 2.4048 the first zero of the Bessel function Jo(z). (bot- 
tom) Maximum value of the population transferred to the 
second qubit in the time interval and with the same parame- 
ters as the figure on top and the following dephasing rates (in 
GHz): 7d cp h = 0.0 (Wue),7dc P h = 0.001 (green) ,j dcp h = 0.005 
(red),7d cp h = 0.01 (turquoise), -y dcph = 0.05 (cyan),-y dcph = 
0.1 (yellow) 



resonance condition is fulfilled and a peak in the popula- 
tion transfer appears induced entirely by the terms in H 2 . 

As a result, the presence of a correction to the canoni- 
cal Hamiltonian Hi does not hinder the manifestation of 
localization effects. By appropriate tuning of the exter- 
nal driving we can select resonance conditions that lead 
to inhibition of transport and provide a fingerprint of the 
underlying coherent evolution. In fig. [4] we illustrate the 
effect of an inhomogeneous distribution of the tunneling 




0.2 0.4 0.6 0.8 
Visibility 



1.0 



FIG. 6: Quantum coherence in a N=2 chain. In this 
graph we plot the magnitude C = max 5^ ^ abs^pij(i)^ as 
a function of the visibility of the transport pattern (population 
fringes in Fig. 5). We have considered 27 values of the pure 
dephasing rate in a logarithmic scale from 7d cp h = 0.1GHz 
(C* m ax — 0) to 7dc P h = 0.0GHz (C max = 0.5). For convenience 
we have fixed 7di ss = 0. 



amplitudes of the qubits within the chain. Not surpris- 
ingly, the presence of this sort of disorder in the array 
leads to a quick loss of contrast. However, the fact that 
the tunneling amplitude of superconducting flux qubits 
is actually tunable can provide mechanism to overcome 
this difficulty by minimizing or even suppressing the local 
disorder in the chain [30]. On the other hand, real im- 
plementations of interacting superconducting flux qubits 
typically result in interaction strengths differences of the 
order of few percents between different nearest neighbors. 
It is important to stress at this point that this differences 
have little effect in our previous discussion since they do 
not affect the resonance conditions in eqs.(15 16), that 



are the key to the hopping inhibition effect presented in 
this work. 

In real implementations of superconducting flux 
qubits we should also expect deviations from a purely 
dissipative model of noise. This motivates the study of 
the effect of pure dephasing processes in our system. 
To this end we will introduce in the master equation 
a Lindblad term of the form given in eq.(22|, with 
dephasing rate 7d ph- To single out the effect of pure 
dephasing, we will take a vanishingly small dissipation 
rate, that is 7di SS = 0. 

In fig. [5] we have studied how the transference of 
population in a chain is affected by pure dephasing 
noise terms. In fig. [BJmp) we plot the dynamics of the 
population for different values of the dephasing rate, 
7dc P h = {0.0, 0.001, 0.01, 0.1}GHz. It is worth remarking 
that the inclusion of these terms in the master equation 
leads to an asymptotic stationary state that is max- 
imally mixed [31 j . However, the transient dynamics 
still provides valuable information. It is remarkable 
to notice that even though the dephasing rapidly 
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scalability of the procedure is illustrated in fig. [7] for an 
array of N — 6 qubits, a system size that is currently 
untractable with tomographic schemes. Notice in this 
graph that, as before, the perfect hopping inhibition is 
only achieved at the zeros of the corresponding Bessel 
functions. The fact that the hopping seems to vanish in 
an extended interval around this point is due to the fi- 
nite size of the time interval considered to perform the 
measure. 



CONCLUSIONS 



FIG. 7: Hopping inhibition in a N = 6 chain. Maximum 
value of the population transferred from the first qubit to the 
last one in a TV = 6 chain within a time interval of 2.8 [is. 
The frequency of the external field is uj — 1.3GHz and its 
magnitude is fixed so that E^/u = 1.2, far from the clos- 
est zero of the Jq(z) Bessel function. The parameters of the 
qubits are Wi = 10GHz and J = lOMHz.For this computation 
we have considered dissipation noise only, with the following 
rates (in GHz): 7 diss = 0.0001 (blue), -y disB = 0.0005 (green), 
Tdiss = 0.001 (red), 7di ss = 0.005 (turquoise, slightly visible). 



destroys the coherent oscillations of the population, 
still the renormalization of the hopping rate survives 
for surprisingly large values of 7d C ph- This fact can be 
clearly appreciated in the panel with E ac /uj = z Q i where 
the hopping is strongly inhibited and still noticeable 
for dephasing rates as high as 7d C ph = 100 MHz. In 
fig. j5jc?oi/;n) we have plotted the same information in 
a more compact format. In this graph the effect of 
dephasing can be clearly seen to reduce the visibility 
of the population transfer oscillations (as a function 
of E ac /uj ). Nonetheless, as stated above, the hopping 
modulation is quite visible also in this graph even for 
large values of the dephasing rate. 

The qualitative relation between the properties of the 
transport dynamics and the coherence in the system 
is illustrated in fig. [6] In this figure we have repre- 
sented the coherence of the chain, quantified by the 
sum of the off diagonal elements of the density matrix 
C = max jjti a bs(pi,j(t)), as a function of the visi- 
bility of the population fringes depicted in fig. [5] For 
the sake of clarity we have used a purely dephasing noise 
since it uniquely affects the coherence terms of the den- 
sity matrix. The main conclusions are however unaffected 
by including dissipation terms. We can see in this graph 
that the degree of quantum coherence of the system can 
be properly quantified by the transport schemes proposed 
in this work. For small values of the dephasing rate, the 
relation between coherence and visibility is indeed essen- 
tially linear. As a result, population measurements alone 
would allow, in the presence of a tunable driving, to de- 
tect signatures of quantum coherence in the system. The 



To summarize, we have analyzed the persistence of 
localization effects beyond exact tight binding Hamiltoni- 
ans and beyond a closed system description. Introducing 
an a.c. interaction term of the form of eq.([2]), we have 
seen that the original ZZ coupling, which provides the 
natural model for the actual coupling in superconduct- 
ing architectures, can be effectively mapped into the 
Hamiltonian eq.(14| where the excitation-preserving 



and non-preserving terms are affected by two different 
renormalized couplings g and g' . We have seen that these 
effective couplings are determined by certain resonance 
conditions (for the coupling g there is always a resonance 
for nf = 0, for g' the stronger condition 2uj + n'ui = 
is required) and, given that the resonance conditions are 
fulfilled, by the arguments of the Bessel functions that 
define these couplings. In conclusion, we have proposed 
a method that allows to tune independently two kinds of 
interaction of very different nature but with the common 
feature of leading to localization phenomena. These are 
shown to be useful for witnessing the coherent behaviour 
of coupled qubit arrays. As a proof of principle, using 
typical parameter regimes in chains of superconducting 
flux qubits, we have shown that transport inhibition 
can be qualitatively linked to the degree of coherence 
in the system. This type of experiments, which involve 
population measurements only, can provide a benchmark 
for quantum behaviour in systems whose complexity 
makes them unsuitcd for detailed tomographic analysis, 
ranging from arrays of self-assembled quantum dots [35] 
to coupled nanomagnets [55] . 
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